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increases. The possibly strong dependence of (qq) on Nf is related to the observed 
Zweig rule violation in the scalar channel. Presuming the existence of a chiral phase 
transition for not too large Nf, we discuss the phenomenological possibilities of sepa- 
rately determining the two-flavour and three-flavour condensates and the quark mass 
/\! ratio r = 2m s /(m u + m^). The issue is closely related to the interpretation of new 

forthcoming precise inr data at low energy. 
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1. Introduction 

The purpose of this paper is to reconsider our understanding of the mechanism of 
spontaneous chiral symmetry breaking (S^SB) in QCD as reflected by the (possibly 
small) size of qq condensates and to interpret its relation to the Zweig rule (ZR) 
violation in the scalar channel and to the quark mass ratio r = 2m s /(m u + m,j). New 
precise experimental results which are expected to come soon [1, 2, 3] will merely 
concern low-energy tttt scattering which is governed by the chiral dynamics of u- and 
rf-quarks with the s-quark playing at most the role of a "sea-side spectator" . From 
this point of view the existence of a relationship between low-energy tttt observables 
and the quark mass ratio r = m s /m (m = m u = rud) is by no means obvious. A 
related question concerns determinations of SU(2) x SU(2) low-energy constants from 
SU(3) x SU(3) observables: for example, the low energy constants Z 4 and / 3 , which 
are essential to assess the prediction of standard chiral perturbation theory (S^PT) 
[4] for tttt s-wave scattering lengths [5], are usually inferred using the experimental 
values of Fx/F n and of the K- and 77-meson masses respectively. 

The theoretical counterpart of these and similar phenomenological questions con- 
cerns the dependence of various order parameters of S^SB on the number Nf of 
massless quark flavours. This ^/-dependence is an effect of light quark loops and it 
is generally expected to be rather weak: it is suppressed in the large- N c limit and 
it violates the Zweig rule which are both considered as a good approximation to the 
real QCD (with a possible exception of the anomalous _+ channel). On the other 
hand, more recent investigations suggest that ^/-dependent light-quark loop effects 
could sometimes be rather important. First, Nature does not seem to respect the 
large- N c predictions in the scalar channel which is not dominated by ideally mixed qq 



nonets, as required by the Zweig rule [6]. Next, some recent lattice simulations with 
dynamical fermions [7, 8] observe a rather strong iV/-dependence of S^SB signals. 
Finally, a new method of estimating the variation of the chiral condensate between 
Nf = 2 and Nf = 3 from the data has recently been proposed [9]. Using as input 
experimental informations on ZR violation in the scalar channel at low and medium 
energies, a large variation of (qq) has been found. 

For sufficiently large Nf/N c , the existence of chiral phase transitions is generally 
expected on the basis of the behaviour of the perturbative QCD /5-function [10]. 
Various approaches have been proposed to study the transitions arising when Nf 
increases. Investigations about the QCD conformal window, where the theory is 
asymptotically free in the ultraviolet, but is governed by a non-trivial fixed point in 
the infrared, suggest a restoration of chiral symmetry for Nf ~ 10 at iV c = 3 [11]. 
The analysis of gap equations in the same conformal window puts a slightly different 
bound, with a critical Nf above 12 [12]. On the other hand, the instanton liquid 
model indicates that instantons do not significantly contribute to Chiral Symmetry 
Breaking for Nf > 6 [13]. 

We would like to stress that particular properties of vector-like gauge theories 
such as QCD allow a new type of non-perturbative interpretation of S^SB and of 
its iV/- dependence, suggesting a natural link between a suppression of (qq) and an 
enhancement of ZR violation in the scalar channel for not too large Nf/N c . No 
reference is made to concepts and methods such as gap equation, mean field ap- 
proximation, or the critical effective coupling, which in confining non-Abelian gauge 
theories lack a truly non-perturbative gauge-invariant definition. Instead, QCD may 
be formulated in Euclidean space-time and the mechanism of S^SB can be related 
to the dynamics of the lowest modes of the hermitean Dirac operator [14, 15, 16] 

H[G\=>y li (d lt + iG li )=Hl[G\, (1.1) 

averaged over gauge field configurations G^(x). In particular, the role of increasing 
Nf naturally appears as a consequence of the paramagnetic effect of light quark loops: 
it tends to suppress the order parameters that are dominated by the small eigenvalues 
of H . Whether this suppression is strong enough to imply a phase transition for not 
too large Nf is a dynamical question which at present we are not able to answer 
analytically. At least one can arrive at a theoretically coherent framework suggesting 
a semi-quantitative understanding of various aspects of S^SB and of their possible 
interconnection. 

2. ^-dependence of infrared-dominated order parameters 

The theory will be considered in a four- dimensional Euclidean box L x L x L x L 
with (anti)periodic boundary conditions modulo a gauge transformation (i.e. on a 



torus). We are interested in gauge invariant correlation functions which break the 
chiral symmetry of the vacuum, i.e. in vacuum expectations of operators which do 
not contain the singlet representation of the chiral symmetry group. All external 
momenta are set to and all quarks are first taken to be massive. At the end we 
consider the limit in which the first Nf lightest quarks become massless keeping the 
remaining quark masses fixed. One randomly chooses a gauge field configuration 
Gn(x) and performs the Grassmann integral over quark fields. The result can be 
formally expressed in terms of eigenvalues and eigenfunctions of the hermitean Dirac 
operator (1.1): 

H[G](f> n = A n [G]0 n , jdx (/&(x)Mx) = S nk- (2.1) 

Since H[G] anticommutes with 75 = 75, the spectrum is symmetric around 0. In 
addition, for gauge field configurations with non-vanishing winding number v G Z, 
the spectrum contains \v\ topological zero modes. Positive eigenvalues will be ranged 
in the ascending order and numerated by a positive integer, defining A_ n = — A n . 
Notice that the eigenvalues A n and the wave functions <f) n (x) are entirely given by 
the gauge field configuration G^x); in particular, they are independent of the quark 
flavour and of the quark masses 1 . The dependence on quark flavour and masses arises 
from the propagator 

S J (x,y\G)=Y: m _ an[G] (2.2) 

and from the fermionic determinant det(M — iH) = TTj A.(rrij\G), where the single 
flavour determinant can be expressed as 

A(m\G)=m^U m2 tt l[ f l - ^ 

n>0 m +UJ n 

The G-independent numbers uj n , which essentially coincide with the free eigenvalues, 
provide a convenient overall normalization of (2.3) and will be specified shortly. An 
integral over fermion fields of any product of bilinear quark currents can be expressed 
in terms of the propagator (2.2) and of the determinant (2.3). The next step then 
consists in taking an average over all gluon configurations. Before we comment on 
this last point, let us concentrate on the simplest example of an order parameter of 
S X SB. 

The chiral condensate (uu), where u stands for the lightest quark, will be con- 
sidered in the SU(N f ) x SU(iV / ) chiral limit 

TOi = m 2 = . . . = mjv,. = m — > 0, rrij 7^ 0, j > Nf (2-4) 



lr The wave functions <j> n live in the spin x colour space and they transform as the fundamental 
representation of SU(iV c ) and as a 4-dimcnsional 0(4) spinor respectively. 



and it will be denoted as — E(iV/). One has 

1 r 1 ttl 

E(N f ) = lim — < / cfeTr 5(x, x|G) > JV/ = lim — < ]T 2 , x2 »*/ • ( 2 - 5 ) 
L 4 j J L 4 „ m^ + A„ ' 

Hereafter, lim denotes the SU(iV/) x SU(iVj) limit (2.4) preceeded by the large volume 
limit. The symbol <^.^$>N f represents the normalized (<C 1 ^>N f — 1) average over 
gauge-field configurations weighted by the fermionic determinant, 

< T > iV/ = Z" 1 /rf/ilG] r A N f(m\G) [] A(m i |G)exp{-5'[G]}, (2.6) 

j>N f 

where S[G] stands for the Yang-Mills action. Since every gauge- field configuration 
Gp(x) can be globally characterized by the corresponding set of Dirac eigenvalues A 
and eigenvectors 0, the functional integral (2.6) may be viewed as an average over 
all possible Dirac spectra. The probability distribution of Dirac eigenvalues should 
be, in principle, calculable from the theory itself. In practice, it requires a non- 
perturbative regularization and renormalization of the gluonic average (2.6) which 
(as in perturbation theory) may depend on the observable I\ Whilst this problem 
is hard to solve in general, in the particular case of chiral order parameters such as 
(2.5), the formal structure of Eq. (2.6) suggests some possibly interesting properties 
of S%SB even before an analytic solution becomes available. 

The fact that (uu) is an order parameter of S^SB is reflected by the vanishing of 
the m — > limit of Eq. (2.5) taken at finite volume. In order to get a non-trivial result, 
the large volume limit has to be taken first and the spectrum of the Dirac operator has 
to become sufficiently dense around the origin. Actually, the average distribution of 
the smallest Dirac eigenvalues is all what matters: if in Eq. (2.5) one cuts the infrared 
end of the spectrum and sticks to |A n | > A, the result would be zero, no matter how 
small A is. For the same reason, the ultraviolet divergences of the sum yj n in Eq. (2.5) 
become irrelevant in the chiral limit. It turns out that only eigenvalues which in the 
average accumulate like 1/L 4 contribute to the chiral condensate (qq) [14, 15, 16]. 
A similar discussion applies to other order parameters of S^SB, in particular to the 
coupling F n of Goldstone bosons to the axial current, whose square can be expressed 
as the two-point correlator of left-handed and right-handed currents (L^R V ) at zero 
momentum transfer. Denoting by F 2 (Nf) the SU(iV/) x SU(Nf) chiral limit (2.4) of 
F%, one has [17] 

F\N S] = lim - « g -j-^-j-^ J fa » N „ (2.7) 



where lim has the same meaning as in Eq. (2.5) and 

1 r 

Jkn= ^ 2 J dx ^l^h^nix) 



(2.8) 



Due to the Goldstone theorem, F 2 (Nf) ^ is both sufficient and necessary for 
the chiral symmetry SU(iV/) x SU(iV/) to be spontaneously broken. Again, F 2 (Nf) 
merely receives contributions from the lowest Dirac eigenvalues but it is less infrared 
sensitive than {uu): the eigenvalues behaving in the average as 1/L 2 could now be 
sufficient to produce a non-zero value of (2.7), since there are two factors m/{m 2 + \ 2 ) 
for a single inverse power of volume [17]. 

Eq. (2.6) suggests that for m — > the effect of the fermionic determinant and the 
iV/-dependence will be stronger for observables T which are dominated by the lowest 
Dirac eigenvalues. This observation follows from the rigorously proven inequality [15] 

l/d 

\\ n [G]\<C—=u n , (2.9) 

where d is the space-time dimension (d — 4 in our case) and C is a constant inde- 
pendent of the gauge field configuration G^(x), of the integer n and of the volume 
V = L d . (In general, C depends on the shape of the space time manifold, once the 
volume has been fixed.) The existence of the uniform upper bound (2.9) reflects the 
paramagnetic response of the Euclidean Dirac spectrum to an external gauge field 
[18]. It allows to split the single flavour determinant (2.3) into infrared and ultravi- 
olet parts [19]: one chooses a cutoff A and defines an integer K such that ujk = A. 
The determinant is then written 

A(m\G) = m lul A m (m\G)A VY (m\G), (2.10) 

where Air involves the first K non-zero eigenvalues and is bounded by 1 as a conse- 
quence of the inequality (2.9), 

airho = n m2 t+ k[ ? ] < 1- (2- 11 ) 

fc=l m +UJ k 

One may expect that in the case of chiral order parameters such as (uu) (2.5) or F 2 
(2.7) which are entirely dominated by the infrared extremity of the Dirac spectrum, 
the Air part of the determinant will describe the bulk of the effect of light quark 
loops in Eq. (2.6). This effect should be paramagnetic, 

£(iV> + 1) < £(iV», F 2 (N f + 1) < F 2 (N f ), (2.12) 

indicating the suppression of chiral order parameters with increasing Nf. How strong 
is this suppression depends on how sensitive is the observable T in Eq. (2.6) to the 
smallest Dirac eigenvalues for large volumes and small quark masses. For this reason 
one may expect a stronger suppression in the case of (qq) than for F n . On the other 
hand, for observables that are not dominated by the infrared extremity of the Dirac 
spectrum, the sensitivity to the determinant and to Nf can remain marginal, as 
expected in the large- N c limit. 



3. Why the scalar channel does not obey large- N c predictions 

We now turn to the connection between flavour dependence of order parameters of 
S^SB and the observed violation of the Zweig rule in the scalar channel. As before, 
we consider first Nf light flavors of common mass m — ► and denote by s the 
(Nf + l)-th quark, whose mass m s is non-zero, but still considered as light compared 
to the scale of the theory. For Nf = 2, this corresponds to the situation in real QCD. 
S(iVj) is a function of m s , and its derivative reads 



d 
dm, 



■£(AM = lim / dx(uu(x)ss(0)) c = U z (m s ) 

m^O J 

=^M?^)(?^H (3i; 



where the notations are as before and the superscript c stands for the connected part. 
Since S(iV/) — > E(7V/ + 1) for m s — > 0, one can write 

X(N f )=X(N f + l) + dfiU z (fi) = J:(N f + l)+m s Z^(m s ) + 0(m 2 s \ogm s ). (3.2) 

«/ u 

In general, and for N f ~ 2 - 3, the difference A(iV/) = S(iV7) - E(7V / + 1) is 
expected to be negligible compared to S(iVj) for two independent reasons: first, it is 
chirally suppressed due to the smallness of m s relative to the condensate S(iV/ + 1), 
at least provided the latter is of the standard "normal" size. Second, the connected 
correlator (3.1) of scalar quark densities of different flavours is suppressed in the 
large- N c limit relative to E(iVj). An important iVj-dependence of the condensate 
would imply that both these arguments fail. We argue that this should naturally be 
expected close to the critical point n cr i t (N c ) at which E(iV/) vanishes. Suppose that 
for a given value of N c (e.g. iV c = 3), we have Nf + 1 < n cr it(A^ c ), so that S(iV/ + 1) 
is still non-zero but already small. Then (for the actual value of m s ) the condensate 
term need not dominate the expansion (3.2) in powers of m s , not because the chiral 
expansion breaks down but due to the suppression of E(iVj + 1). The generalized 
chiral perturbation theory (G;\TT) is precisely designed to cope with such a situation 
[20, 21, 22]. The suppression of the condensate means that near the critical point, 
the average density of states [14] 

-^ « P (A) > >iV/ , p(A) = ^ S(X - Xn[G]), (3.3) 

drops for A ~ m. Indeed, Eq. (2.5) can be rewritten as 

E(iV» = 21imy o —— 2 -^ p {mu)^ Nf . (3.4) 



It is then natural to expect that the proximity of a phase transition will further 
manifest itself by an increase of fluctuations of the density of states and/or by an 



enhancement of the density- density correlation L <Cp(A)p(A') ^>% for A ~ A' ~ m. 
This quantity determines the variation of the condensate, see Eqs. (3.1) and (3.2): 

9 „/, r x . ,. f°° du dv 1 , . , . . , n . 

_E(iV / ) = 4b m y o —— - <<R(mu)R(m , v)>>%r (3.5) 

Far away from the critical point, the correlation of small Dirac eigenvalues should not 
be important, as predicted by the large- N c limit. Since this limit leads to a suppres- 
sion of quark loops, any infrared-dominated order parameter becomes independent 
of Nf. Therefore, the large-iV c limit prevents £ from vanishing, since the critical 
number of flavours n cr it(N c ) moves away to infinity as N c — > oo. This asymptotic 
behaviour is also supported by perturbative calculations, in which Nf and N c usu- 
ally arise through their ratio. Large- N c expansion is therefore expected to converge 
slowly for Nf fixed just below the critical point n cr it(N c ), and to yield irreparably 
false results for Nf above it. This argument merely concerns the vacuum channel 
++ . The variation of any chiral order parameter between Nf and Nf + 1 is given 
by a correlation function which violates the Zweig rule precisely in that channel, cf 
Eq. (3.2): a strong variation would imply the existence of J PC = ++ , SLV(iVf)- 
singlet states strongly coupled both to the first Nf light quarks and to the scalar 
density ss of the (Nf + l)-th quark. The proximity of a phase transition could then 
explain in a natural way why the spectrum of ++ states is not dominated by ideally 
mixed scalar mesons, presenting significant discrepancies with large- N c predictions 2 . 
In the following we concentrate on the actual case Nf = 2 — 3, having in mind the 
possibility that (for N c = 3) the real world might already be close to the critical point. 
It has recently been pointed out [9] that the correlation function (3.1) satisfies a well 
convergent sum rule which allows a phenomenological estimate of the iVj-dependence 
of the condensate between Nf = 2 and Nf = 3, 

Tl z (m s ) = - I™ -^ a(t), (3.6) 

n Jo t 

where the spectral function a(p 2 ) (defined in Minkowski space-time) collects ZR 
violating contributions which couple both to uu and ss 

<P 2 ) = iE( 2 ^ (4) (p-^)(0|^|n)(n|s S |0). (3.7) 

n 

Since we take m u = rrid = m and isospin symmetry cannot be spontaneously bro- 
ken [24], only isoscalar states \n) contribute in Eq. (3.7) and (0\uu\n) = (0\dd\n). 
Under the plausible assumption that the dominant contribution comes from two- 
particle states (\irir),\KK), . . .), the spectral function (3.7) can be reconstructed 
from the corresponding multi-channel T-matrix, which contains the experimental 

2 Connections with scenarios invoking trace anomaly and light dilatons [23] remain to be seen. 
We thank H. Leutwylcr for bringing our attention to this question. 



information on the size of Zweig rule violation in the ++ channel, such as the 
effect of the /o(980) resonance. Using further \PT to normalize the solution of 
the multi-channel Muskhelishvili-Omnes equation, one can obtain the corresponding 
form factors {0\qq\n). Within the standard version of %PT, a strong iVj-dependence 
of the condensate has been found [9]: S(3)/S(2) = 1 — 0.54 ± 0.27. We shall further 
comment on this result shortly. 

4. xPT considerations 

The case of a large Zweig-rule violation leading to a substantial difference between 
the Nf = 2 and Nf = 3 condensates has never been fully included into the %PT 
analysis before. This leads us to carefully reconsider the G;\RT relation between 
the actual size of the condensate(s), the quark mass ratio r = m s /m, and some low 
energy observables. What matters is the renormalization group invariant product 
mE in physical units of F%M%, i.e. the Gell-Mann-Oakes-Renner ratio(s) [25] 

7T 7T 

The standard chiral expansion presumes that both X(2) and X(3) are close to 1 
which is, a priori, hardly compatible with an important flavour dependence of the 
condensate. On the other hand, in G%PT, the GOR ratio (4.1) can be significantly 
below 1, leaving enough space for a large iVf -dependence. It is worth stressing that 
an important ZR violation would not affect the quantitative relation between X(2) 
and low energy nn observables provided it is systematically based on SU(2) x SU(2) 
G%PT. The reason is not a "bad convergence" of the expansion in the strange quark 
mass. The ZR violation effects do not manifest themselves through a specific low- 
energy constants in the Nf = 2 effective Lagrangian, whereas in the three-flavour 
£ c ff they show up as extra (ZR violating) low energy constants Lq(h) and L^fj,) 
which are a priori unknown and have never been determined experimentally 3 . Hence, 
X(2) remains accessible to experiment via low-energy im phases [21, 22], azimuthal 
asymmetries in the decay r — > 3 r nv T [26], or in the reaction 7 + 7 — ► 3ti [27]. The 
question remains how this information is related to the size of X(3), to the ZR 
violation in the ++ channel and to the quark mass ratio r = m s /m. 

A partial answer to these questions can be obtained from a new look at old 
SU(3) x SU(3) G%PT expansions of kaon and pion masses (see Ref. [20] and App. A 
of Ref. [22]), rewritten as 



FlMl = 2m 



X(3) + (m s + 2m)Z* s (m s ) + 4m 2 A eS + F^M^ (4.2) 



3 The reason of this difference resides in group theory: while the ZR violating correlator 
(uu(x)dd(O)} does not break the chiral symmetry SU(2)xSU(2), the two-point function (uu(x)ss(O)) 
is an order parameter for both Nf — 2 and Nf — 3 chiral symmetry. 



f k m k = (m + m s ) E(3) + (m s + 2m)Z^ e (m s 



+ (m + m s ) 2 A cS + Fl5M 



L K- 



(4.3) 



These formulae collect in a scale independent manner all contributions linear and 
quadratic in quark masses. They are useful to the extent that 5Mp <C Mj>, which 
is certainly a much weaker requirement than the assumption which is at the basis 
of S%PT. The constants Z^ s and A c r are related to the quark-mass independent 



constants of the 0(p 2 ) G%PT Lagrangian, Zq and Aq, renormalized at a scale [i: 
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(4.5) 



where we use the notations B = E(3)/F 2 (3), Ml v = lim m ^ ^l, v - Both expres- 
sions (4.4) and (4.5) are scale independent. The connection with the S%PT 0(p 4 ) 
constants is F 2 (3)Z,f(/i) = 16B 2 L 6 (fi) and F 2 (3)A (n) = 16B 2 L 8 (fi) respectively. 
The constant Z^(m s ), independent of m u = m d = m, is the same as in Eq. (3.2), 
taken for Nf = 2. It is convenient to split the scale independent remainders 5M 2 K 
into two scale independent parts, 5M 2 = 5(\)M 2 + 5{2)M 2 : 



S(D M * 
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(4.7) 



Substituting Eqs. (4.6) and (4.7) in Eqs. (4.2) and (4.3), one recovers the full 0(p 4 ) 
standard expansion [4]. In this case, 5(2) M 2 consists of 0(p & ) (two- loop) and higher 
standard contributions [28]. In the G%PT reading, <5( 2 )M 2 are 0(m^ ulirk ) as well, 
but now, they consist of the tree £(0,3) component of C c g, of the remaining scale 
independent part of the one-loop 0(p 4 ) contributions not included in Eqs. (4.6) and 
(4.7), and of higher order terms. As a result we expect 5M 2 K /M 2 K to be at most 
3-4 per cent in the whole range < X(3) < 1 (this statement will be made more 
quantitative in the final result). The control of the accuracy, independently of the 
size of X(3), is considerably simplified expanding the product FpMp rather than 
Mp and Fp separately. It avoids uncertainties related to the low energy constant £ 



(1/5 ) an d to its ZR violating counterpart £ (I/4). Further advantages of this way of 
ordering the expansion of Goldstone boson masses will appear shortly. 

A simple algebra allows one to infer from Eqs. (4.2) and (4.3) a relation between 
X(3), the ZR violating constant Z^ s and the quark mass ratio r = m s /m: 

x(3) + 2m(m. + 2m) z& = x _ _ (r) + ^^ (4g) 

where SX(3) is a simple combination 4 of 5M% K and 

7T 7T 

This information can now be combined with the general Eq. (3.2) (considered here for 
Nf = 2). The latter can be recovered considering the limit S(2) = \im m ^o(F^M^) /2m 
of Eq. (4.2) keeping m s fixed: 

X{2) = X(3) + jg^Zi + SX(2). (4.10) 

Eliminating Z^ s from Eqs. (4.8) and (4.10), one arrives at a simple relation between 
the Nf = 2 and Nf = 3 GOR ratios (4.1) and the quark mass ratio r 

X(2) = [1 - e(r)]-J- + -4oX(3) + A. (4.11) 

r + 2 r + 2 

Before we show that the remainder A is small and well under control, the simple 
meaning of Eq. (4.11) should be stressed. The three-flavour GOR ratio should be 
in the interval < -X"(3) < X{2) because of the vacuum stability and of the para- 
magnetic inequality (2.12). Since, furthermore, r > f\ = 2F K M K /F 7C M n — 1 ~ 8, 
the second term on the r.h.s. of Eq. (4.11) represents a small correction all over the 
interval of r: the quark mass ratio m s /m is merely given by the two-flavour GOR 
ratio X(2) which is more easily accessible experimentally, whereas X(3) represents 
the fine tuning of the above relation. It might be, for instance, conceivable that 
X(2) would be close to 1, the quark- mass ratio r close to its standard value (or even 
larger) and yet -X"(3) ~ 0, implying a very steep decrease of X(Nf) to the critical 
point and a huge amount of ZR violation. This is a possibility which has never been 
considered before. 

We finally discuss the uncertainty in Eqs. (4.8), (4.10) and (4.11). One has 



SX(3) 



e r 



r — 1 



5M 2 K r + \5Ml 
M\ r - 1 Ml ' 



(4.12) 



whereas SX{2) can be read off from the expansion of F%M% inspecting linear terms 
in the limit m — > with m s fixed (including the whole 0(p 4 ) G%PT order, see 

4 Hereaftcr, all manipulations with Eqs. (4.2) and (4.3) are algebraically exact, making no use of 
expansions in 6M% K or in quark masses. 
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Eq. (5.46) in [22]). As a result, SX(2) = 5 {2 )M%/M% + . . ., where the dots stand for 
contributions for which the dimensional estimate gives 0[m 2 m s / {M^K.h)\ ~ ±0.4/r 2 . 
Hence, this uncertainty is at the per mille level. A is a simple linear combination 
A = 5X (2) + 5X (3) -r / (r + 2) and it may be represented as A = Ai±e. Ai represents 
the contributions of chiral logs to (4.12) arising from <5(i)M 2 K . The contribution of 
5(2)M 2 to A almost exactly cancels between 8X{3) and 8X(2) and the remaining 
uncertainty e is dominated by S^M^. The dimensional estimate gives e = 0.025 for 
r = 10, e = 0.006 for r = 20 and even much smaller values for higher r. In Fig. 1, 
the correlation between X(2) and r is plotted, including A^ Eq. (4.10) is considered 
either with a maximal ZR violation [A(3) = 0] or no violation at all [A(3) = A(2)]. 
We check the fine tuning role of X(3) in the significant correlation between r and 
X(2). If the two-flavour GOR ratio is close to 1, the quark mass ratio is not very 
accurately determined, but it is much more restricted for smaller X(2). 
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Figure 1: Correlation between the Gell-Mann-Oakes-Renner ratio for two massless 
flavours X(2) and the quark mass ratio r = m s /m. The uncertainty due to e is not 
taken into account (see text). The upper curve respects strictly the Zweig rule, which is 
maximally violated on the lower curve. 



Further insight can be obtained combining the general perturbative expressions 
displayed above with Moussallam's sum rule (3.6). Differentiating Eq. (3.2) with 
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respect to m s and using Eq. (4.4) yields the identity 

2mm 



X(2)-X(3) 



<<s 



FlMl 

IX 7T 

where Xp = m s • d(\og Mp)/dm SJ and 



n z (m s ) + ^- 2 (x K + lx v 



+ AX, (4.13) 



MIFl \ dm J m->o 2m v ; 



*■ , 



The identity (4.13) is a variant of Eq. (4.10), in which the ZR violating constant Z^ s 
has not been eliminated but reexpressed using the sum rule (3.6) taken at m = 0. AX 
differs from <5X(2) by insertion of (l—m s d/dm s ). They are both expected of the same 
order of magnitude and of opposite sign. In particular, AX receives contributions 
starting at the two-loop order in S%PT. Following the same dimensional estimate 
as in the discussion of the uncertainty in Eq. (4.11), one expects |AX| < 0.05 for 
r ~ 10, further decreasing for larger r. 

Eq. (4.13) provides a general framework for the discussion of the variation of 
the condensate X(Nf) between Nf = 2 and Nf = 3 in terms of the sum rule (3.6), 
independently of any prejudice about the size of X(3). It is in particular valid in 
S%PT, in which the deviation of X(3) from 1 is considered as a small perturbation. 
In the latter case, A^ = A,, = 1 at the leading order, and 2mm s I?Q can be replaced by 
M^[M\ — M%/2). In this way one recovers the S^PT-based analysis of Ref. [9]. The 
second term on the right hand side of Eq. (4.13) then takes the value 0.21, whereas 
the evaluation of the sum rule (3.6), as discussed in Ref. [9], corresponds to the final 
result 0.38 < X(2) — X(3) < 0.73, which is compatible with the conclusion expressed 
in Ref. [9] in terms of the ratio X(3)/X(2). Notice that Uz(m s ) > for not too large 
m s , because the two-point function (3.1) exhibits a positive logarithmic increase as 
m s — ► 0. Consequently, in S%PT, the bound X(2) — X(3) > 0.21 must hold up to 
two-loop corrections. A too large difference X(2) — X(3) could hardly be reconciled 
with the premises of S;\PT which require both GOR ratios to be reasonably close 
to 1. Consequently, a new analysis of the sum rule (4.13) within G%PT would be 
highly desirable. 

5. Summary and concluding remarks 

1. Order parameters of SB^S which are dominated by the infrared extremity of the 
Euclidean Dirac spectrum (E = —(uu), F%, ...) could be sensitive to the para- 
magnetic effect of light quark loops : in the chiral limit, the fermionic determinant 
reduces the statistical weight of the lowest Dirac eigenvalues and gradually suppresses 
these order parameters as the number of massless flavours Nf increases. This might 
result into a rich chiral phase structure as a function of Nf and N c . 
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2. For Nf approaching the critical point n cr i t (N c ), the average density of low 
Dirac eigenvalues should drop and its fluctuations should increase. This would natu- 
rally lead to a significant reduction of (qq), i.e. to a rapid decrease of the Gell-Mann- 
Oakes-Renner (GOR) ratio X(Nf) (4.1), and to an enhancement of the Zweig rule 
(ZR) violation in the scalar channel as compared to the large- N c predictions. 

3. In the real world, where m u ~ m^ <C m s <C A# ~ 1 GeV, X{2) and X(3) are 
of a direct phenomenological interest (1 > X(2) > X(3) > 0). If the critical point 
(for N c =3) is far away from Nf = 3, it is natural to expect X(3) ~ X(2) ~ 1 and a 
less important ZR violation in the ++ channel. If, on the other hand, the real world 
is close to a phase transition, X(Nf) should quickly fall towards the critical point, 
leading to a large ZR violating difference X(2) — X(3). This would force a value of 
X(3) significantly below 1, leaving open the question whether X(2) already feels the 
influence of the critical point or still remains close to 1. 

4. X(2) can be extracted from precise low-energy 7T7r-scattering experiments, 
independently of the ZR violation and of the size of X(3). Furthermore, X(2) is 
closely related to the quark mass ratio r = 2m s /(m u + m^), and this relation is 
only marginally affected by X(3). On the other hand, even if X(2) ~ 1, and S%PT 
were a reliable expansion scheme in the two-flavour sector, its accurate predictions 
for s-wave scattering lengths [4, 5] would be obstructed by important ZR violation 
already at the one-loop level: sofar, a reliable determination of the SU(2) x SU(2) 
low-energy 0(p A ) constants Z3 and U from independent experimental data requires 
(among other things) the knowledge of the SU(3) x SU(3) ZR violating constants L e 
and L 4 . The determination of the two-flavour GOR ratio X(2) remains the central 
goal of ongoing irir scattering experiments [1, 2, 3] and of related proposals [26, 27]. 
If X{2) turned out to be close to 1, these experiments could be interpreted as a first 
measurement of the low-energy constants I3 and I4. 

5. Whatever the experimental output for X(2) will be, additional information 
will be necessary to pin down X(3) and to settle the theoretical issue of a nearby phase 
transition. The sum rule analysis of Ref. [9] could be extended, but a more direct 
access to the ZR violation in the ++ channel and to the three-flavour condensate 
would be desirable despite its difficulty 

6. The theoretical question of what happens in the phase in which we likely 
do not live [Nf > n CT i t (N c )\, is at present hard to answer unambiguously. The first 
chiral transition should merely affect observables that are particularly sensitive to 
the lowest modes of the Dirac operator. The p-meson mass, the string tension and 
the characteristics of confinement in general need not be affected at all. Among chiral 
order parameters, (qq) exhibits the strongest infrared sensitivity and is expected to 
vanish first. Whether this implies a complete or partial restoration of chiral symmetry 
[29] is a question that involves the hard problem of non-perturbative regularization 
and renormalization of the chiral symmetry-breaking sector of QCD. Within the cut- 
off dependent bare theory, the vanishing of (qq) implies F% = 0, and consequently 
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the full restoration of chiral symmetry [30] . The validity of this argument in the full 
theory remains to be clarified. The answer likely resides in a non-perturbative study 
of the fixed points of renormalization group flows. 
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